Doppler wind LIDAR systems data processing and applications : an overview towards developing the new generation of wind remote-sensing sensors for off-shore wind farms

This Ph.D. thesis addresses remote sensing of the atmosphere by means of lidar and S-band clear-air weather radar, and related data signal processing. Active remote sensing by means of these instruments offers unprecedented capabilities of spatial and temporal resolutions for vertical atmospheric profiling and the retrieval of key optical and physical atmospheric products in an increasing environmental regulatory framework. The first goal is this Ph.D. concerns the estimation of error bounds in the inversion of the profile of the atmospheric backscatter coefficient from elastic lidar signals (i.e., without wavelength shift in reception when interacting with atmospheric scatterers) by means of the two-component inversion algorithm (the so-called Klett-Fernald-Sasano’s algorithm). This objective departs from previous works at the Remote Sensing Lab. (RSLab) of the Universitat Politècnica de Catalunya (UPC) and derives first-order error-propagated bounds (approximate) and total-increment bounds (exact). As distinctive feature in the state of the art, the error bounds merge into a single body both systematic (i.e., user-calibration inputs) and random error sources (finite signal-to-noise ratio, SNR) yielding an explicit mathematical form. The second goal, central to this Ph.D., tackles retrieval of the Atmospheric Boundary Layer Height (ABLH) from elastic lidar and S-band Frequency-Modulated Continuous-Wave (FMCW) radar observations by using adaptive techniques based on the Extended Kalman Filter (EKF). The filter is based on morphological modelling of the Mixing-Layer-to-Free-Troposphere transition and continuous estimation of the noise covariance information. In the lidar-EKF realization the proposed technique is shown to outperform classic ABLH estimators such as those based on derivative techniques, thresholded decision, or the variance centroid method. The EKF formulation is applied to both ceilometer and UPC lidar records in high- and low-SNR scenes. The lidar-EKF approach is re-formulated and successfully extended to S-band radar scenes (Bragg’s scattering) in presence of interferent noise sources (Rayleigh scattering from e.g., insects and birds). In this context, the FMCW feature enables the range-resolved capability. EKF-lidar and EKF-radar ABLH estimates are cross-examined from field campaign results. Finally, the third goal deals with exploitation of the existing UPC lidar station: In a first introductory part, a modified algorithm for enhancing the dynamic range of elastic lidar channels by “gluing” analog and photon-counting data records is formulated. In a second part, two case examples (including application of the gluing algorithm) are presented to illustrate the capabilities of the UPC lidar in networked atmospheric observation of two recent volcano eruption events as part of the EARLINET (European Aerosol Research Lidar Network). The latter is part of GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems) framework.La tesis doctoral aborda la teledetecció atmosfèrica amb tècniques lidar i radar (banda S) i llur tractament del senyal. La teledetecció activa amb aquests instruments ofereix resolucions espacials i temporals sense precedents en la perfilometria vertical de l’atmosfera i recuperació de productes de dades òptics i físics atmosfèrics en un marc de creixent regulació mediambiental. El primer objectiu d’aquesta tesi concerneix l’estimació de cotes d’error en la inversió del perfil del coeficient de retrodispersió atmosfèrica a partir de senyals lidar de tipus elàstic (és a dir, sense desplaçament de la longitud d’ona en recepció al interactuar amb els dispersors atmosfèrics) mitjançant l’algorisme d’inversió de dues components de Klett-Fernald-Sasano. Aquest objectiu parteix de treballs previs en el Remote Sensing Lab. (RSLab) de la Universitat Politècnica de Catalunya (UPC) i permet obtenir cotes de primer ordre (aproximades) basades en propagació d’errors i cotes (exactes) basades en el increment total de l’error. Característica diferencial en front l’estat de l’art és l’assimilació d’errors sistemàtics (per exemple, entrades de cal·libració d’usuari) i aleatoris (relació senyal-soroll, SNR, finita) en forma matemàtica explícita. El segon objectiu, central de la tesis, aborda l’estimació de l’altura de la capa límit atmosfèrica (ABLH) a partir de senyal lidar elàstics i d’observacions radar en banda S (ona continua amb modulació en freqüència, FMCW) utilitzant tècniques adaptatives basades en filtrat estès de Kalman (EKF). El filtre es basa en modelat morfològic de la transició atmosfèrica entre la capa de mescla i la troposfera lliure i en l’estimació continua de la informació de covariança del soroll. En el prototipus lidar-EKF la tècnica proposada millora clarament les tècniques clàssiques d’estimació de la ABLH como són les basades en mètodes derivatius, decisió de llindar, o el mètode de la variança-centroide. La formulació EKF s’aplica tant a mesures procedents de ceilòmetres lidar como de la pròpia estació lidar UPC en escenes d’alta i baixa SNR. Addicionalment, l’enfoc lidar-EKF es reformula i s’estén amb èxit a escenes radar en banda S (dispersió Bragg) en presència de fonts de soroll interferent (dispersió Rayleigh de, per exemple, insectes i ocells). En aquest context, la característica FMCW permet la capacitat de resolució en distància. L’estimació de la ABLH amb els prototipus lidar-EKF i radar-EKF s’interrompés en campanyes de mesura. Finalment, el tercer objectiu atén a l’explotació de l’estació lidar UPC existent: En una primera part introductòria, es formula un algorisme modificat de “gluing” per a la millora del marge dinàmic de canals lidar elàstics mitjançant combinació (o "enganxat") de senyals lidar adquirits analògicament i amb foto-comptatge. En una segona part, es presenten dos exemples (incloent l’aplicació de l’algorisme de “gluing”) que il·lustren les capacitats del lidar de la UPC en l’observació atmosfèrica de dos recents erupcions volcàniques des de la xarxa d’observació EARLINET (European Aerosol Research Lidar Network). Aquesta última és part de GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems).Postprint (published version

Doppler wind LIDAR systems data processing and applications : an overview towards developing the new generation of wind remote-sensing sensors for off-shore wind farms

This Ph.D. thesis addresses remote sensing of the atmosphere by means of lidar and S-band clear-air weather radar, and related data signal processing. Active remote sensing by means of these instruments offers unprecedented capabilities of spatial and temporal resolutions for vertical atmospheric profiling and the retrieval of key optical and physical atmospheric products in an increasing environmental regulatory framework. The first goal is this Ph.D. concerns the estimation of error bounds in the inversion of the profile of the atmospheric backscatter coefficient from elastic lidar signals (i.e., without wavelength shift in reception when interacting with atmospheric scatterers) by means of the two-component inversion algorithm (the so-called Klett-Fernald-Sasano’s algorithm). This objective departs from previous works at the Remote Sensing Lab. (RSLab) of the Universitat Politècnica de Catalunya (UPC) and derives first-order error-propagated bounds (approximate) and total-increment bounds (exact). As distinctive feature in the state of the art, the error bounds merge into a single body both systematic (i.e., user-calibration inputs) and random error sources (finite signal-to-noise ratio, SNR) yielding an explicit mathematical form. The second goal, central to this Ph.D., tackles retrieval of the Atmospheric Boundary Layer Height (ABLH) from elastic lidar and S-band Frequency-Modulated Continuous-Wave (FMCW) radar observations by using adaptive techniques based on the Extended Kalman Filter (EKF). The filter is based on morphological modelling of the Mixing-Layer-to-Free-Troposphere transition and continuous estimation of the noise covariance information. In the lidar-EKF realization the proposed technique is shown to outperform classic ABLH estimators such as those based on derivative techniques, thresholded decision, or the variance centroid method. The EKF formulation is applied to both ceilometer and UPC lidar records in high- and low-SNR scenes. The lidar-EKF approach is re-formulated and successfully extended to S-band radar scenes (Bragg’s scattering) in presence of interferent noise sources (Rayleigh scattering from e.g., insects and birds). In this context, the FMCW feature enables the range-resolved capability. EKF-lidar and EKF-radar ABLH estimates are cross-examined from field campaign results. Finally, the third goal deals with exploitation of the existing UPC lidar station: In a first introductory part, a modified algorithm for enhancing the dynamic range of elastic lidar channels by “gluing” analog and photon-counting data records is formulated. In a second part, two case examples (including application of the gluing algorithm) are presented to illustrate the capabilities of the UPC lidar in networked atmospheric observation of two recent volcano eruption events as part of the EARLINET (European Aerosol Research Lidar Network). The latter is part of GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems) framework.La tesis doctoral aborda la teledetecció atmosfèrica amb tècniques lidar i radar (banda S) i llur tractament del senyal. La teledetecció activa amb aquests instruments ofereix resolucions espacials i temporals sense precedents en la perfilometria vertical de l’atmosfera i recuperació de productes de dades òptics i físics atmosfèrics en un marc de creixent regulació mediambiental. El primer objectiu d’aquesta tesi concerneix l’estimació de cotes d’error en la inversió del perfil del coeficient de retrodispersió atmosfèrica a partir de senyals lidar de tipus elàstic (és a dir, sense desplaçament de la longitud d’ona en recepció al interactuar amb els dispersors atmosfèrics) mitjançant l’algorisme d’inversió de dues components de Klett-Fernald-Sasano. Aquest objectiu parteix de treballs previs en el Remote Sensing Lab. (RSLab) de la Universitat Politècnica de Catalunya (UPC) i permet obtenir cotes de primer ordre (aproximades) basades en propagació d’errors i cotes (exactes) basades en el increment total de l’error. Característica diferencial en front l’estat de l’art és l’assimilació d’errors sistemàtics (per exemple, entrades de cal·libració d’usuari) i aleatoris (relació senyal-soroll, SNR, finita) en forma matemàtica explícita. El segon objectiu, central de la tesis, aborda l’estimació de l’altura de la capa límit atmosfèrica (ABLH) a partir de senyal lidar elàstics i d’observacions radar en banda S (ona continua amb modulació en freqüència, FMCW) utilitzant tècniques adaptatives basades en filtrat estès de Kalman (EKF). El filtre es basa en modelat morfològic de la transició atmosfèrica entre la capa de mescla i la troposfera lliure i en l’estimació continua de la informació de covariança del soroll. En el prototipus lidar-EKF la tècnica proposada millora clarament les tècniques clàssiques d’estimació de la ABLH como són les basades en mètodes derivatius, decisió de llindar, o el mètode de la variança-centroide. La formulació EKF s’aplica tant a mesures procedents de ceilòmetres lidar como de la pròpia estació lidar UPC en escenes d’alta i baixa SNR. Addicionalment, l’enfoc lidar-EKF es reformula i s’estén amb èxit a escenes radar en banda S (dispersió Bragg) en presència de fonts de soroll interferent (dispersió Rayleigh de, per exemple, insectes i ocells). En aquest context, la característica FMCW permet la capacitat de resolució en distància. L’estimació de la ABLH amb els prototipus lidar-EKF i radar-EKF s’interrompés en campanyes de mesura. Finalment, el tercer objectiu atén a l’explotació de l’estació lidar UPC existent: En una primera part introductòria, es formula un algorisme modificat de “gluing” per a la millora del marge dinàmic de canals lidar elàstics mitjançant combinació (o "enganxat") de senyals lidar adquirits analògicament i amb foto-comptatge. En una segona part, es presenten dos exemples (incloent l’aplicació de l’algorisme de “gluing”) que il·lustren les capacitats del lidar de la UPC en l’observació atmosfèrica de dos recents erupcions volcàniques des de la xarxa d’observació EARLINET (European Aerosol Research Lidar Network). Aquesta última és part de GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems)

Vertical Azimuth Display simulator for wind-Doppler lidar error assessment

(c) 2017 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes, creating new collective works for resale or redistribution to servers or lists, or reuse of any copyrighted components of this work in other works.This works presents a simplified Vertical Azimuth Display (VAD) motion simulator for off-shore wind lidars. The simulator is rooted to the case of a conically-scanning lidar (e.g., the Zephyr lidar), where the wind speed vector is retrieved from the Line-of-Sight velocities over one scan period. The methodological part addresses the geometrical foundations of the simulator and how the lidar attitude is assimilated in matrix form. The discussion part considers the case of time-invariant, horizontally-homogeneous wind under two motional cases of the lidar, static and dynamic. Cases examples are parameterized by Horizontal Wind Speed, Wind Direction and tilt amplitude.Postprint (author's final draft

Estimation of the motion-induced horizontal-wind-speed standard deviation in an offshore doppler lidar

This work presents a new methodology to estimate the motion-induced standard deviation and related turbulence intensity on the retrieved horizontal wind speed by means of the velocity-azimuth-display algorithm applied to the conical scanning pattern of a floating Doppler lidar. The method considers a ZephIR™300 continuous-wave focusable Doppler lidar and does not require access to individual line-of-sight radial-wind information along the scanning pattern. The method combines a software-based velocity-azimuth-display and motion simulator and a statistical recursive procedure to estimate the horizontal wind speed standard deviation—as a well as the turbulence intensity—due to floating lidar buoy motion. The motion-induced error is estimated from the simulator’s side by using basic motional parameters, namely, roll/pitch angular amplitude and period of the floating lidar buoy, as well as reference wind speed and direction measurements at the study height. The impact of buoy motion on the retrieved wind speed and related standard deviation is compared against a reference sonic anemometer and a reference fixed lidar over a 60-day period at the IJmuiden test site (the Netherlands). Individual case examples and an analysis of the overall campaign are presented. After the correction, the mean deviation in the horizontal wind speed standard deviation between the reference and the floating lidar was improved by about 70%, from 0.14 m/s (uncorrected) to -0.04 m/s (corrected), which makes evident the goodness of the method. Equivalently, the error on the estimated turbulence intensity (3–20 m/s range) reduced from 38% (uncorrected) to 4% (corrected).Peer ReviewedPostprint (published version

The Use of Augmented Reality to Strengthen Competence in Data Analysis and Problem Solving in Engineering Students at the Universidad del Valle de México

The objective of this research was to analyze the improvement in the data analysis and problem-solving competence of students of industrial and systems engineering (IIS) and mechatronics engineering (IMEC) through the use of this technology and its impact on the results of the undergraduate general examination (EGEL). A training course was held for teachers and students for the design of learning objects (LO), and a questionnaire on the use of AR and the improvement in learning was administered. AR is a technology that has begun to be introduced in different contexts and at different educational levels. The results obtained through the Wilcoxon test and the multiple correspondence analysis (MCA) showed that there were improvements in academic performance with the use of AR and an interest in this tool being used during the academic training process.Laureate EducationRevisión por pare

Doppler wind LIDAR systems data processing and applications : an overview towards developing the new generation of wind remote-sensing sensors for off-shore wind farms

This Ph.D. thesis addresses remote sensing of the atmosphere by means of lidar and S-band clear-air weather radar, and related data signal processing. Active remote sensing by means of these instruments offers unprecedented capabilities of spatial and temporal resolutions for vertical atmospheric profiling and the retrieval of key optical and physical atmospheric products in an increasing environmental regulatory framework. The first goal is this Ph.D. concerns the estimation of error bounds in the inversion of the profile of the atmospheric backscatter coefficient from elastic lidar signals (i.e., without wavelength shift in reception when interacting with atmospheric scatterers) by means of the two-component inversion algorithm (the so-called Klett-Fernald-Sasano’s algorithm). This objective departs from previous works at the Remote Sensing Lab. (RSLab) of the Universitat Politècnica de Catalunya (UPC) and derives first-order error-propagated bounds (approximate) and total-increment bounds (exact). As distinctive feature in the state of the art, the error bounds merge into a single body both systematic (i.e., user-calibration inputs) and random error sources (finite signal-to-noise ratio, SNR) yielding an explicit mathematical form. The second goal, central to this Ph.D., tackles retrieval of the Atmospheric Boundary Layer Height (ABLH) from elastic lidar and S-band Frequency-Modulated Continuous-Wave (FMCW) radar observations by using adaptive techniques based on the Extended Kalman Filter (EKF). The filter is based on morphological modelling of the Mixing-Layer-to-Free-Troposphere transition and continuous estimation of the noise covariance information. In the lidar-EKF realization the proposed technique is shown to outperform classic ABLH estimators such as those based on derivative techniques, thresholded decision, or the variance centroid method. The EKF formulation is applied to both ceilometer and UPC lidar records in high- and low-SNR scenes. The lidar-EKF approach is re-formulated and successfully extended to S-band radar scenes (Bragg’s scattering) in presence of interferent noise sources (Rayleigh scattering from e.g., insects and birds). In this context, the FMCW feature enables the range-resolved capability. EKF-lidar and EKF-radar ABLH estimates are cross-examined from field campaign results. Finally, the third goal deals with exploitation of the existing UPC lidar station: In a first introductory part, a modified algorithm for enhancing the dynamic range of elastic lidar channels by “gluing” analog and photon-counting data records is formulated. In a second part, two case examples (including application of the gluing algorithm) are presented to illustrate the capabilities of the UPC lidar in networked atmospheric observation of two recent volcano eruption events as part of the EARLINET (European Aerosol Research Lidar Network). The latter is part of GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems) framework.La tesis doctoral aborda la teledetecció atmosfèrica amb tècniques lidar i radar (banda S) i llur tractament del senyal. La teledetecció activa amb aquests instruments ofereix resolucions espacials i temporals sense precedents en la perfilometria vertical de l’atmosfera i recuperació de productes de dades òptics i físics atmosfèrics en un marc de creixent regulació mediambiental. El primer objectiu d’aquesta tesi concerneix l’estimació de cotes d’error en la inversió del perfil del coeficient de retrodispersió atmosfèrica a partir de senyals lidar de tipus elàstic (és a dir, sense desplaçament de la longitud d’ona en recepció al interactuar amb els dispersors atmosfèrics) mitjançant l’algorisme d’inversió de dues components de Klett-Fernald-Sasano. Aquest objectiu parteix de treballs previs en el Remote Sensing Lab. (RSLab) de la Universitat Politècnica de Catalunya (UPC) i permet obtenir cotes de primer ordre (aproximades) basades en propagació d’errors i cotes (exactes) basades en el increment total de l’error. Característica diferencial en front l’estat de l’art és l’assimilació d’errors sistemàtics (per exemple, entrades de cal·libració d’usuari) i aleatoris (relació senyal-soroll, SNR, finita) en forma matemàtica explícita. El segon objectiu, central de la tesis, aborda l’estimació de l’altura de la capa límit atmosfèrica (ABLH) a partir de senyal lidar elàstics i d’observacions radar en banda S (ona continua amb modulació en freqüència, FMCW) utilitzant tècniques adaptatives basades en filtrat estès de Kalman (EKF). El filtre es basa en modelat morfològic de la transició atmosfèrica entre la capa de mescla i la troposfera lliure i en l’estimació continua de la informació de covariança del soroll. En el prototipus lidar-EKF la tècnica proposada millora clarament les tècniques clàssiques d’estimació de la ABLH como són les basades en mètodes derivatius, decisió de llindar, o el mètode de la variança-centroide. La formulació EKF s’aplica tant a mesures procedents de ceilòmetres lidar como de la pròpia estació lidar UPC en escenes d’alta i baixa SNR. Addicionalment, l’enfoc lidar-EKF es reformula i s’estén amb èxit a escenes radar en banda S (dispersió Bragg) en presència de fonts de soroll interferent (dispersió Rayleigh de, per exemple, insectes i ocells). En aquest context, la característica FMCW permet la capacitat de resolució en distància. L’estimació de la ABLH amb els prototipus lidar-EKF i radar-EKF s’interrompés en campanyes de mesura. Finalment, el tercer objectiu atén a l’explotació de l’estació lidar UPC existent: En una primera part introductòria, es formula un algorisme modificat de “gluing” per a la millora del marge dinàmic de canals lidar elàstics mitjançant combinació (o "enganxat") de senyals lidar adquirits analògicament i amb foto-comptatge. En una segona part, es presenten dos exemples (incloent l’aplicació de l’algorisme de “gluing”) que il·lustren les capacitats del lidar de la UPC en l’observació atmosfèrica de dos recents erupcions volcàniques des de la xarxa d’observació EARLINET (European Aerosol Research Lidar Network). Aquesta última és part de GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems)

Performance evaluation of a floating lidar buoy in nearshore conditions

This work provides a signal-processing and statistical-error analysis methodology to assess key performance indicators for a floating Doppler wind lidar. The study introduces the raw-to-clean data processing chain, error assessment indicators and key performance indicators, as well as two filtering methods at post-processing level to alleviate the impact of angular motion and spatial variability of the wind flow on the performance indicators. Towards this aim, the study mainly revisits horizontal wind speed (HWS) and turbulence intensity measurements with a floating ZephIR 300 lidar buoy during a 38 day nearshore test campaign in Pont del Petroli (Barcelona). Typical day cases along with overall statistics for the whole campaign are discussed to illustrate the methodology and processing tools developed.Peer Reviewe

Vertical Azimuth Display simulator for wind-Doppler lidar error assessment

(c) 2017 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes, creating new collective works for resale or redistribution to servers or lists, or reuse of any copyrighted components of this work in other works.This works presents a simplified Vertical Azimuth Display (VAD) motion simulator for off-shore wind lidars. The simulator is rooted to the case of a conically-scanning lidar (e.g., the Zephyr lidar), where the wind speed vector is retrieved from the Line-of-Sight velocities over one scan period. The methodological part addresses the geometrical foundations of the simulator and how the lidar attitude is assimilated in matrix form. The discussion part considers the case of time-invariant, horizontally-homogeneous wind under two motional cases of the lidar, static and dynamic. Cases examples are parameterized by Horizontal Wind Speed, Wind Direction and tilt amplitude

Estimation of wave period from pitch and roll of a lidar buoy

This work proposes a new wave-period estimation (L-dB) method based on the power-spectral-density (PSD) estimation of pitch and roll motional time series of a Doppler wind lidar buoy under the assumption of small angles (±22 deg) and slow yaw drifts (1 min), and the neglection of translational motion. We revisit the buoy’s simplified two-degrees-of-freedom (2-DoF) motional model and formulate the PSD associated with the eigenaxis tilt of the lidar buoy, which was modelled as a complex-number random process. From this, we present the L-dB method, which estimates the wave period as the average wavelength associated to the cutoff frequency span at which the spectral components drop off L decibels from the peak level. In the framework of the IJmuiden campaign (North Sea, 29 March–17 June 2015), the L-dB method is compared in reference to most common oceanographic wave-period estimation methods by using a TriaxysTM buoy. Parametric analysis showed good agreement (correlation coefficient, ρ = 0.86, root-mean-square error (RMSE) = 0.46 s, and mean difference, MD = 0.02 s) between the proposed L-dB method and the oceanographic zero-crossing method when the threshold L was set at 8 dB.This research was funded by the Spanish Government and EU Regional Development Funds, ARS project PGC2018-094132-B-I00 and ACTRIS-IMP project GA-871115. The European Institute of Innovation and Technology (EIT), KIC InnoEnergy project NEPTUNE (call FP7) supported the measurement campaigns. The Generalitat de Catalunya—AGAUR funded doctoral grant 2020 FISDU 00455 by A. Salcedo-Bosch. CommSensLab-UPC is an Excellence Unit (MDM-2016-0600) funded by the Agencia Estatal de Investigación, Spain.Peer ReviewedPostprint (published version

Floating Doppler Wind Lidar measurement of wind turbulence: a cluster analysis

The standard deviation of the Horizontal Wind Speed as a proxy of wind turbulence is used to compare the apparent wind turbulence measured by an off-shore floating Doppler lidar to the one measured by a fixed lidar on a metmast. We use statistical analysis based on clustering the horizontal wind speed measured by the floating lidar as well as buoy angular amplitude and period under the approximation of harmonic motion. Three scenarios with different wave and wind conditions are discussed from the IJmuiden's test campaign (North Sea.).This work was supported via Spanish Government–European Regional Development Funds project PGC2018-094132-B-I00 and EU H2020 ACTRIS-2 (GA 654109). The European Institute of Innovation and Technology (EIT), KIC InnoEnergy project NEPTUNE (Offshore Metocean Data Measuring Equipment and Wind, Wave and Current Analysis and Forecasting Software, call FP7) supported measurements campaigns. CommSensLab is a María-de-Maeztu Unit of Excellence funded by the Agencia Estatal de Investigacion (Spanish National Science Foundation).Peer ReviewedPostprint (author's final draft